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SECTION  1 
INTRODUCTION 


Adaptive  equalization  is  one  of  several  mitigation  schemes  that 
can  be  designed  into  HF  communication  systems  to  combat  the  effects  of 
frequency  selective  fading.  Other  possible  schemes  include  maximum  like¬ 
lihood  sequence  estimation  and  various  spread  spectrum  techniques.  The 
restriction  on  the  scope  of  this  work  to  adaptive  equalization  should  not 
be  interpreted  as  a  disparaging  reflection  upon  the  potential  benefits  of 
the  other  possibilities,  but  merely  represents  the  limitations  caused  by 
time  and  resource  constraints. 

The  use  of  adaptive  equalization  does  not  increase  the  bandwidth 
of  the  transmitted  signal.  Also,  for  the  narrowband  modulation  type  of 
interest  here,  PSK,  long,  deep  fades  can  occur  which  span  the  modulation 
bandwidth.  Thus  one  would  envision  an  adaptive  equalizer  to  be  used  in 
the  context  of  HF  transceivers  having  the  capability  to  scan  the  HF  band 
and  to  select  and  maintain  a  good  carrier  frequency  automatically. 

The  particular  type  of  equalizer  addressed  in  this  report  is 
called  a  decision  feedback  equalizer  (DFE),  so  named  because  decisions  on 
the  value  of  previously  detected  channel  symbols  are  fed  back  to  the 
equalizer  to  be  incorporated  into  its  internal  state  vector. 

Section  2  of  this  report  is  meant  to  be  a  complete  introduction 
to  adaptive  decision  feedback  equalization.  By  collecting  the  pertinent 
contents  of  several  journal  articles  into  this  section  the  reader  is  freed 
from  sorting  through  several  sets  of  notation.  Also  all  derivations  in 
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Section  2  are  generalized  for  fractionally-spaced  sampling,  coherent 
demodulation  and  time-variant  channels,  while  the  journal  articles  tend  to 
address  special  cases.  Of  the  numerous  articles  on  adaptive  equalization 
to  be  found  in  the  open  literature,  three  in  particular  are  widely  refer¬ 
enced,  historically  important,  and  pertinent  to  this  report:  Austin's 
derivation  of  the  decision  feedback  equalizer  from  decision-theoretic 
principles  (Reference  1);  Gersho's  treatment  of  the  stochastic  gradient 
algorithm  in  the  context  of  linear  equalization  (Reference  2)  which  is 
easily  extended  to  decision  feedback  equalization;  and  Godard's  applica¬ 
tion  of  Kalman  filtering  principles  to  equalizer  adaptation  (Reference 
3).  In  addition  to  these.  References  4,  5,  and  6  develop  the  fast  Kalman 
adaptation  algorithm.  The  derivation  of  the  fast  Kalman  algorithm  in 
Section  2  closely  follows  the  treatment  of  transversal  filter  equalizer 
adaptation  found  in  Reference  6,  but  extends  the  derivation  to  decision 
feedback  equalizer  adaptation  using  the  mathematical  tools  provided  by 
Reference  4. 

The  results  obtained  during  this  investigation  were  generated 
with  a  Monte  Carlo  computer  simulation  of  an  HF  link.  The  computer  simu¬ 
lation  is  described  in  Section  3,  with  primary  emphasis  given  to  the 
channel  and  receiver  models. 

Section  4  discusses  the  results  of  a  simulation-aided  analysis 
of  HF  PSK  modems.  First,  the  performance  of  conventional  and 
equalizer-assisted  modems  are  compared  for  a  simple  time-invariant  HF 
channel  model.  Then  an  equalizer-assisted  modem  is  analyzed  for  a  simple 
time-varying  frequency-selective  HF  channel  model. 
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SECTION  2 

ADAPTIVE  DECISION  FEEDBACK  EQUALIZATION 


A  decision  feedback  equalizer  (DFE)  is  a  digital  signal  proces¬ 
sor  that  can  enhance  the  performance  of  several  types  of  digital  modems 
operating  under  frequency  selective  channel  conditions.  This  discussion 
is  generalized  to  include  T/N  sampling,  phase  coherent  demodulation  and 
random  channel  variability.  T/N  sampling  means  the  receiver  samples  the 
analog  channel  output  every  T/N  seconds  where  T  is  the  minimum  signaling 
interval  of  the  modulated  waveform.  Phase  coherent  demodulation  means  the 
demodulator  uses  phase  as  well  as  amplitude  information,  which  implies  the 
DFE  must  process  complex  numerical  data.  Random  channel  variability  means 
the  channel  transfer  function  is  time-variant  and  is  unknown  a  priori  to 
the  receiver,  which  implies  the  DFE  must  be  able  to  continuously  adapt  its 
internal  parameters  to  match  the  channel. 

Two  adaptation  techniques  commonly  discussed  in  the  technical 
literature  are  the  fast  Kalman  and  the  stochastic  gradient  algorithms. 

Both  of  the  algorithms  have  larger  implementation  requirements  than  the 
DFE  itself.  The  fast  Kalman  algorithm  requires  more  hardware  and/or 
digital  processing  than  that  required  to  implement  the  stochastic  gradient 
algorithm.  However,  the  fast  Kalman  algorithm  converges  faster  than  the 
stochastic  gradient  algorithm. 

We  shall  first  discuss  a  baseband  channel  representation  that  is 
appropriate  for  several  digital  modulation  schemes,  and  then  describe  the 
T/N  fractionally-spaced  DFE.  The  remainder  of  this  report  formulates 

i 

several  adaptation  techniques  starting  with  the  stochastic  gradient 
algorithm.  Then  various  least  square  adaptation  algorithms  shall  be 
derived,  culminating  with  the  fast  Kalman  algorithm. 


2.1 


BASEBAND  CHANNEL  MODEL 


It  is  convenient  to  model  the  modulation,  channel  propagation 
and  demodulation  processes  with  a  time-discrete  complex  baseband  (or 
analytic  signal)  representation.  Such  a  representation  of  the  modulated 
waveform  is  the  sequence  of  complex  samples  {d|<}.  The  transmitted  wave¬ 
form  represented  by  the  discrete  sequence  { d ^ }  is  of  the  form 


V  [ Rea  1  (d.  )  cos(m  t)  +  Imag(d,  )  sin(u)  t)]  rect  (t/T  -  k) 

li  K  C  K  C  S 

k=0 

where  u>c  is  the  carrier  frequency,  Ts  is  the  sampling  interval,  and  the 
rect  function  is  defined  as 

I  1  -1/2  <  X  £  1/2 

rect  (X)  =  ' 

(  0  otherwise  . 

Several  quadrature  amplitude  modulation  (QAM)  and  M-ary  phase 
shift  keying  (PSK)  modulation  types  can  be  represented  with  this  discrete 
baseband  formulation  by  appropriately  selecting  the  alphabet  for  d^ .  For 
example,  limiting  d^  to  the  set  { 1 , - 1 }  would  correspond  to  BPSk  modu¬ 
lation.  Also,  T/N  spaced  sampling  is  represented  by  only  changing  d^ 
every  N  samples,  in  which  case  T^  =  T/N.  We  shall  use  the  con.v.ntion  that 
the  receiver  generates  a  detected  symbol  when  k  is  an  integer  multiple  of 
N. 


The  time-discrete  analytic  signal  representation  of  the  channel 
at  the  kth  sample  time  is  the  set  of  complex  samples  {h™},  the  sampled- 
dat.a  response  of  the  channel  to  an  isolated  unit  sample.  More  specifical¬ 
ly,  suppose  that  dn  is  one  and  d|<  is  zero  for  all  k  not  equal  to  n.  Then 
the  discrete  channel  output  (ignoring  noise)  at  the  kth  sample  time  is 


O 

V 


k  -n-M 

h^  .  The  integer  M  corresponds  to  the  propagation  delay  and  is  used 
here  as  a  convenience  so  that  we  can  require  h™  equal  to  zero  for  m  <  0, 
since  the  channel  is  causal  and  hence  M  must  be  non-negative.  It  is 
further  assumed  that  the  channel  has  finite  memory,  and  we  introduce  L, 
called  the  channel  memory  duration,  which  is  the  smallest  integer  such 
that  h1^.1  is  zero  for  all  m  >  L.  For  our  purposes  the  values  h1^  in  the 
range  0  <  m  <  L  completely  specify  the  intersymbol  interference  character¬ 
istics  of  the  channel  and  are  used  to  form  the  vector  h^ . 

4  [hk  hk  hk  •••t,k)- 


The  underbar  denotes  a  column  vector  and  superscript  T  denotes 
vector  transpose.  The  transpose  form  is  used  for  typographical  simplicity 
because  it  allows  the  column  vectors  to  be  presented  in  a  horizontal 
orientation.  Next  we  define  a  vector  of  the  most  recent  L+l  channel 
symbol s . 

|T 


4 


fdk  dk-l  dk  -2"  ’  -l  1 


k-U 


(1) 


The  kth  received  sample  constituting  the  output  of  the  discrete 
channel  ,  denoted  uk ,  is  the  convolution  of  the  unit  sample  response  with  a 
sequence  of  inputs,  corrupted  by  additive  noise. 


hk  dk-m  *  \ 

m  =  U 


(2) 


m,.  is  the  complex  white  gaussian  noise  component  received  during 
the  kt*1  sampling  interval.  <\  block  diagram  of  the  tapped  delay  line 
implementation  of  the  baseband  channel  model  is  shown  in  Figure  1. 


1)  A  linear  forward  predictor  of  the  next  N+l  samples  to  enter 
xj<;  is  derived.  This  forward  predictor  uses  the  same  Kalman 
gain  vector  of  an  LS  algorithm  adapting  a  DFE. 

2)  A  linear  backward  predictor  of  the  last  N+l  samples  to  leave 
xj<  is  derived.  This  backward  predictor  also  uses  ^ . 

3)  LS  adaptation  is  applied  to  an  extended  state  DFE  to  develop 
the  mathematical  structure  between  the  forward  and  backward 
predictors. 

4)  It  is  shown  that  the  forward  and  backward  predictors  can 
work  in  tandem  to  efficiently  calculate  ^'recursively. 

2.6.1  Forward  Prediction 


Define  a  vector  ^  composed  of  the  N+l  elements  entering  the 
,tate  vector  at  the  kth  sample  time  when  k/N  is  an  integer. 


,T  a  r 

A  =  >uk  Vi 


"k-N+l  dk-D-N 


The  linear  prediction  of  _§j<+n  is  F|<_n  xj<  where  F|<  is  an  (N+l) 
x  (K i +1+K2 )  rectangular  matrix  of  forward  prediction  coefficients.  The 
forward  prediction  error  vector  is  defined  as 


f  =  f  _  p  V 

-4c  -VN  r k-N  * 


(40) 


A  comparison  of  this  equation  with  Equation  (21)  shows  a 
similarity  between  this  forward  predictor  and  the  DFE.  To  calculate  the 
optimal  value  of  Fk,  the  scalar  error  LS  cost  function  of  Equation  (8)  is 


The  sequence  of  operations  performed  by  this  algorithm  after  the 
k^  sample  time  (when  k/N  is  an  integer)  is 


^k  =  Tk-N  h 

(34) 

[33,34] 

3k  =  V(X  +  V 

(35) 

[30,34, 

35] 

Tk  =  “  tTk-N  "  h 

(36) 

[21] 

£k  =  dk-D  '  VN  \ 

(37) 

[23] 

3k  =  3k-N  +  3]<  ek  • 

(38) 

The  implementation  complexity  of  this  algorithm  can  be  dominated 
by  the  first  step  which  requires  (K1+I+K2)2  multiplications.  The  follow¬ 
ing  algorithm  is  equivalent  to  this  one,  but  is  computationally  superior 
because  it  manipulates  smaller  matrices. 

2.6  THE  FAST  KALMAN  ADAPTATION  ALGORITHM 

None  of  the  adaptation  algorithms  discussed  thus  far  take  into 
consideration  any  deterministic  relationship  between  xj<  and  xj<-N*  From 
the  definition  of  xy  in  Equation  (5),  all  but  N+l  of  the  elements  in  21k -N 
also  reside  in  The  fast  Kalman  algorithm  attains  its  computational 

efficiency  by  exploiting  this  fact. 

The  derivation  of  the  fast  Kalman  Algorithm  will  proceed  as 

fol lows : 
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Since  X  +  1S  a  scalar  and  is  invertible. 


A-1  =  —  A'1 

k  X  k-N 


Ak-N  ^k  *-k  Ak-N 
X  +  4  Ak-N  4 


Define  Tk  =  A^  and  notice  that  Tk  like  Ak  is  hermitian  since 


Tk  Ak  =  (AkTk)T  =  1  =  (AkV*  =  Tk  AK  =  TK  A 


T  T  * 

and  postmultiplication  by  the  inverse  of  Ak  implies  that  Tk  equals  Tk< 

Substituting  Tk  into  Equation  (29)  yields  the  following  formula. 


T  -  I  T  -  Tk-N  4  h  Tk-N 

k  "  x  k"N  "  x'+>  T 

4  'k-N  ^k 


Now  consider  g,  . 


[19]  %  '  \  ik  *  Tk  *k 


[30,31] 


T  x 
k-N  ~k 


x  *  h  Tk-«4 


T  x 
k-N  -k 


A  +  *k  Tk-N  *k 
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2.5  THE  KALMAN  ADAPTATION  ALGORITHM 

First  we  shall  derive  a  recursive  calculation  of  from 
which  a  recursive  formula  for  cy<  will  easily  follow.  This  algorithm 
bypasses  the  necessity  of  matix  inversion,  and  as  a  result  is  more  effi¬ 
cient  than  the  one  discussed  immediately  above. 

It  is  easy  to  verify  the  following  equation. 


Ak-N  h  X  ^k  =  Ak-N  h  X  +  h  Ak-N 

r  *  -1  t1  * 

[*  +  4Ak_N£*J  ^ 


*  .-1 


Ak-N  ^k  A  h  =  Ak-N  h.  fA  +  h  Ak-N  4  + 

*-l  ..  ,-l  *  a-1  ..  ,  .. 


Ak-N  h  X  h  Ak-N  fA  +  Ak-N  ^k 


i  ■  [i +  \1N  ^  **]  u  -  +  h  \-N  -xj,)'1  x**] 


k-N  "k-N  -A-k ;  ^k 


Postmultiply  by  Ak_N  and  premultiply  by  Ak_N 


1  "  tAk-N  "  ±k  X  ^k ^  ^Ak-N 


Ak-N  ^k  (A  +  *k  Ak-N  ^  ^k  A k -N -I  {'27>l 


[13,27]  I  =  I  A.  [A'1,  -  A"*  x.  (A  +  x*  A^  xj’1  x*  A'1 


k  L  k-N  k-N  ^k 


^k  k-N  H< ;  ^k  k-N j 


mm 


VC:  v*v 


The  recursive  calculation  of  c^  is  now  derived. 


[10,14] 

Ak  4 

=  X  Ak-N  4-N  +  4  dk-D 

(15) 

[13,15] 

Ak  4 

=  tAk  ’  4  4^  4-N  +  4  dk-D 

(16) 

★  -k 

[16] 

Ak  4 

=  \  4-n  +  4  ^dk-D  "  4-n  4^ 

(17) 

Premultiply  by  A"1 

[U]  ^  ^  (dR_Q  -  c^_N  XjJ  (18) 

where  ^  =  A^1  ^  •  (19) 

The  sequence  of  operations  of  this  recursive  algorithm  after  the 


kth  sample  time  is 

[13]  Ak  =  A  Ar_n  +  2^  (20) 

ek  =  dk-D  "  4-N  4 

[19]  ^  =  A”1  (22) 

[18,21]  ^  _N  +  ^  er  .  (23) 


Notice  that  the  calculation  of  A^  and  cy  is  greatly  simplified 
and  that  the  calculation  of  the  quantity  yj<  has  been  bypassed  altogether. 
However,  A^  must  be  inverted  to  calculate  ^  and  the  dimension  of  A|<  is 
Kj+l+Kj,  often  a  large  number.  The  calculation  of  can  also  be  done 
recursively,  and  algorithms  that  do  so  are  discussed  in  the  next  two 
sections . 
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The  quantity  on  the  right  is  a  column  vector  and  the  quantity 
premultiplying  c^  on  the  left  is  a  positive-definite  hermitian  square 
matrix.  Some  of  the  forthcoming  derivations  are  tedious.  To  simplify  the 
presentations,  the  numbers  of  preceding  equations  justifying  a  new 
expression  are  listed  in  brackets  to  the  left,  as  was  done  above.  Equa¬ 
tion  (9)  can  be  restated  as 


[9] 

\  4s  4 

(10) 

where 

[9,10] 

A  =  T  l"-k'N  X  x* 

Ak  niQ  ^nN  *fiN 

(11) 

[9,10] 

n-k/N  * 

^k  =  ^  X  ^fiN  dnN-D 

n=u 

(12) 

The  recursive  nature  of  Equations  (11)  and  (12) 

is  apparent 

can  be 

written  as 

[11] 

flk  "  1  flk-N  *  2*  h 

(13) 

[12] 

Xk  ■  *  ik-N  +  4  dk-D 

(14) 

A  new  coefficient  vector  is  calculated  after  each  equalizer 
operation.  Thus  at  the  kth  sample  time  the  quantities  A^,  vj<_n  and  cj<_n 
are  available  so  that  A^ ,  v]<  and  c|<  may  be  calculated  recursively. 


2.4 


LEAST  SQUARES  ADAPTATION  AND  THE  WIENER -HOPF  EQUATIONS 


The  stochastic  gradient  algorithm  discussed  above  was  based  upon 
the  least  mean  square  (LMS)  cost  function.  Now  we  turn  our  attention  to 
algorithms  that  minimize  the  least  squares  (I.S)  cost  function  given  by 


LS  cost  function 


k/N 


I 

n  =  0 


n-  k/N 

A 


dnN-D 


4-x-nN 


(8) 


This  cost  function  is  the  geometrically  weighted  sum  of  the  squared 
magnitude  of  all  equalizer  output  errors  if  the  equalizer  coefficient 
vector  had  always  been  equal  to  cj<.  After  the  kth  equalizer  output  has 
been  generated,  the  least  squares  adaptation  algorithms  generate  the  value 
of  Ck  that  minimizes  the  LS  cost  function.  As  done  previously,  the 
discrete  independent  temporal  index  k  increments  with  each  received  sample 
and  hence  increases  by  N  with  each  equalizer  operation.  All  quantities 
discussed  henceforth  are  defined  only  when  k/N  is  an  integer.  Thus,  the 
previous  value  of  quantities  that  change  only  after  each  equalizer  opera¬ 
tion  are  denoted  with  the  subscript  k-N  instead  of  k-1. 

For  time  invariant  channels  the  selection  of  cj<  should  depend  as 
much  on  previous  data  as  current  data,  in  which  case  the  constant  para¬ 
meter  x  appearing  in  the  cost  function  should  be  unity.  X  is  incorporated 
into  the  LS  cost  function  to  allow  for  time  varying  channels.  A  real 
positive  value  of  X  slightly  smaller  than  unity  causes  received  data  to  be 
less  influential  as  it  becomes  older. 


The  optimal  value  of  c^  is  found  by  setting  the  gradient  of  the 
!.S  cost  function  with  respect  to  cj<  to  zero  which  yields  the  discrete 
Wiener-Mopf  equations. 


$ 


Each  term  in  this  equation  can  be  readily  derived. 


ik-N  ik 


=  2  x*  x*  Cj<_N  =  2 


*  ^dk-D  - k-N  h.  +  dk-D  ^k-N  ^  ^  =  2  ^ 


k-D  4 


Therefore 


5  hi2  ■  -2  (dk-o  -4>4  ■ -2  u 


Then  the  stochastic  gradient  adaptation  rule  becomes: 


£k  ♦  Sk.N  *  6  £k  £k 


To  state  the  operation  of  a  T/N  fractionally-spaced  DFE  using 
the  stochastic  gradient  adaptation  algorithm  we  define  a  square  matrix  and 
a  column  vector  both  of  dimension  K1+I+K2  and  a  rectangular  matrix  having 
dimensions  (K1+I+K2)  x  N 


4  I 


K^l-N  ! 


^  (IC1+2)th  element 


a  =  [0  0  ...  0  1  0  ...  0] 


All  unlabeled  partitions  contain  zero  elements  only.  I  denotes  the  iden¬ 
tity  matrix  with  dimension  specified  by  the  subscript.  The  adaptive 
equalizer  operation  is  given  by  Algorithm  1. 


This  shall  be  done  using  a  numerical  search  technique  based  upon  the 
method  of  steepest  descent  called  the  stochastic  gradient  algorithm.  The 
stochastic  gradient  iteration  rule  is 

c*  «-  ck-N  -  1/2  A  v| ek | 2  .  (6) 

V j  ek J  2  denotes  the  approximation  of  the  gradient  of  |ek|2  with 
respect  to  £k-N  where  the  slight  dependence  of  cj<_n  on  d^-D  and  is 

ignored.  A  is  a  real  scalar  controlling  algorithm  performance,  and  the 
factor  1/2  is  introduced  here  for  later  rotational  convenience.  The  use 
of  the  term  "gradient"  conforms  to  the  nomenclature  of  common  engineering 
practice  and  the  technical  literature.  However,  such  usage  is  mathemati¬ 
cally  imprecise  because  the  gradient  of  any  real-valued  function  with 
respect  to  a  complex  valued  vector  either  does  not  exist  or  is  identically 
zero.  To  eliminate  confusion  the  reader  must  interpret  the  gradient  of  a 
real -valued  function  with  respect  to  a  complex-valued  vector  as  the  sum  of 
1)  the  gradient  of  the  real -valued  function  with  respect  to  the  real  part 
of  the  complex  vector  and  2)  the  product  of  the  square  root  of  -1  and  the 
gradient  of  the  real -valued  function  with  respect  to  the  imaginary  part  of 
the  complex  vector. 

II  2 

£|<  3t  |  points  in  the  direction  of  maximum 
increase  of  | | z .  The  stochastic  gradient  algorithm  attempts  to  minimize 
| | 2  by  perturbing  c^_^  in  the  direction  approximately  opposite  that  of 
this  gradient  after  each  equalizer  operation.  To  derive  V  Je^j2,  first 
expand  |e^j2. 


9|  £k  |  =  9  1  ^dk-D  '  ^k-N  ^dk-D  “  ^k-N  ^ 


9  I dk-D  2  +  9  |^k-N  h\2  + 


"  9  K-D  Ck-N  ^k  +  dk-D  ^k-N  ^  ^ 


+ 


not  have  a  priori  knowledge  of  dk_D-N»  in  which  case  the  decision  symbol 
dk-D-N  is  based  upon  the  equalizer  output  yk-N*  However,  to  simplify  our 
consideration  of  the  adaptation  algorithms  we  shall  assume  that  knowledge 
of  dk-o-N  is  available  to  the  receiver.  The  complex  tap  weights  of  the 
feedforward  section  are  chosen  to  coherently  combine  all  the  signal 
components  of  the  received  data  arising  from  the  transmission  of  dk-D  and 
can  be  viewed  as  a  transversal  filter  qualitatively  matched  to  the 
discrete  channel  unit  sample  response.  This  linear  combination  of  the 
received  data  symbols  will  also  have  signal  components  arising  from 
transmitted  channel  symbols  other  than  dk-D»  half  of  which  were  transmit¬ 
ted  prior  to  and  half  after  dk_D*  Presumably,  those  transmitted  prior  to 
dk-D  have  already  been  detected  correctly.  Hence  the  associated  inter¬ 
symbol  interference  ( I S I )  components  in  the  weighted  sum  of  the  received 
data  are  deterministically  related  to  the  previous  decision  symbols.  The 
coefficients  of  the  feedback  section,  then,  are  selected  to  cancel  these 
components.  For  the  T/N  fractionally-spaced  DFE,  L/N  previously  detected 
symhols  are  contributing  ISI  components  to  the  weighted  sum  of  the 
feedforward  taps,  and  hence  KL  should  be  set  equal  to  L  and  K2  should  be 
set  approximately  equal  to  L/N. 

2.3  STOCHASTIC  GRADIENT  ALGORITHM 

Define  the  error  in  the  kth  equalizer  output  as  ^  =  dk_D  -  yk» 
where  the  availability  of  dk-D  is  assumed.  One  would  like  to  iterate  cj< 
to  eventually  drive  ^k  to  zero,  but  cannot  completely  diminish  due  to 
the  random  noise  component  in  yk-  Failing  that,  it  is  desirable  and 
mathematically  convenient  to  have  the  adaptation  algorithm  try  to  minimize 
the  least  mean  square  (LMS)  cost  function. 

LMS  cost  function  =  E { | ek  J  2 } 


diagram  of  T/N  fractionally-spaced  decision  feedback 


2.2 


THE  T/N  FRACTIONALLY  SP*CED  OFE 


According  to  the  definitions  given  above,  a  particular  transmit¬ 
ted  channel  symbol  dk_|_  influences  the  currently  received  channel  symbol 

uk  and  previously  received  symbols  uk_i,  'Jk-2 .  Uk-L*  equalizer 

ideally  uses  all  of  these  received  symbols  to  assist  the  receiver  in 
generating  the  decision  symbol  dk-L»  the  receiver's  estimate  of  dk-L»  As 
mentioned  earlier,  the  values  of  d^  change  at  most  every  N  samples,  and 
the  receiver  generates  a  decision  symbol  only  when  k/N  is  an  integer.  To 
do  this  the  equalizer  actually  processes  the  currently  received  symbol  and 
the  Ki  previously  received  symbols,  where  Ki  is  chosen  to  be  greater  than 
or  equal  to  L.  A  OFE  processes  K2  past  decision  symbols  also.  After  the 
kth  sample  time  the  decision  dk-D  is  made,  and  the  integer  0  must  he  in 
the  range  L  <  D  Kj  for  proper  equalizer  operation,  yk  denotes  the  OFE 
output  generated  after  the  k.th  sample  and  is  calculated  by 


yk  =  ^<-M  ^ 


(*  denotes  conjugate  transpose)  (3) 


where 


*  A 
= 


r  -Kl  -Ki  +1 

Lck  ck 


r-i  r0  ri 
ck  ck  'k 


<?I 


(4) 


and 


T 

^k 


4  K 


Uk-1  Uk-Ki 


dk-D-N  dk-0-2rr”dk-0-NK2^ 


(5) 


The  mechanization  of  this  calculation  is  shown  in  Figure  2.  xj,.  is  called 
the  state  vector  and  c^.m  the  coefficient  vector.  The  elements  of  c^-N 

are  the  tap  weights  shown  in  the  block  diagram.  For  this  coherent 
baseband  implementation  of  the  equalizer,  these  coefficients  are  complex 
sea  1 ars  . 


The  block  diagram  i no!  :ates  the  DFE  can  be  implemented  as  two 
tapped  delay  lines.  The  one  on  the  left  is  called  the  feedforward  section 
and  the  one  on  the  right,  the  feedback  section.  Usually  the  receiver  does 


discrete  baseband  channel 


forward  prediction  =  l  X  lInN+N-Fk*flN ^inN+N"Fk-*nN 

LS  cost  function  n=0  "  k^N  ^  N  k^N 


Each  column  of  F|<  contains  all  the  coefficients  used  to  predict 
a  single  element  within  The  Wiener-Hopf  equations  yielding  the 

optimal  F|<  are  found  by  setting  to  zero  the  gradient  of  the  trace  of  the 
LS  cost  function  with  respect  to  each  column  of  F^ .  The  result  is 


k/N  .  ... 
(  l  x"-k/N 
n=0 


*  If  -  k(N  xn_k/N 
^oN^nN-^k  ~  ^  X 

n=0 


^-nN  — nN 


[11,42] 


[42,43] 


where 

f  J  kf  xn-k/N 

k  n 
n=0 


*fiN  — nN+N 


X  Vk-N  +  4  ik+N  ‘ 


Repeating  the  manipulations  of  Equations  (15)  through  (19)  using 
F|<_n  and  j]<+N  instead  of  cy-N  and  d|<_D  yields 

f18]  Fk  =  Fk-N  *  9J4+N  '  Fk-N4^  *45* 

[40.45]  F|<  »  Fk_K  *  ^  4  (46) 


Now  we  define  the  minimum  value  of  the  LS  cost  function  given  by 
Equation  (41)  as  E  which  is  found  by  substituting  Equation  (43)  into 
Equation  (41)  and  simplifying. 


,"»*  «£  *  * 
’*  1%  *  **  L  •  C  *  L 


k  ~  L 
K  n=0 


^nN+N  -nN+N 


‘k-N  -^k+N  ^k+N 


A  recursive  calculation  of  E|<  is  now  derived. 

[4M«-  Efc  ■  x  Rk-N  *  4+ N  4+N  -  K-Nt2*4*N>*(Fk-NtSkfk)  <49) 

47 , 4o  J 

C4M6.  Ek  -  x(«k-N-4f>k-N)  *  (I*+Fk-A)(fk+Fk-Nik)* 

-  4  *{!„  %  4*  -  4+n  4*  Fk  (60> 

C-'O  44  Ekf  .  X  ff_N  *  4tN  x*(Fk_N-Fk)  +44**  F*_n  h  4* 

47 , bO J 

-  (vkf*-  ik.Ni**)  4  4*  (61) 

[46,51]  Ekf  .  X  ef_n  +  4+n  4(-44)  +  (4  +  F*_  4  -  vk*  4)  4 


+  ^k+N  h,  h 

(52) 

[52] 

Ek  =  x  eJ-n  +  (Ik+  Fk-N^k“Vk  ^  ^k 

(53) 

[43.46] 

MFk-N+9k-k)  =  Vk 

(54) 

[54] 

h  4  =  Ak*  Vk  '  Fk-N 

(55) 

[19,55] 

f-k  h  h  =  Vk*  2*  -  Fk-N  h 

(56) 

★ 

Since  ^  1S  a  scalar, 

[53,56]  [[  -  ie[+  U-44)  ^  4  •  (57) 

One  more  identity  is  needed  for  the  forward  predictor. 

:4°,56]  4  -  4tN  *  4  4  i)c  •  vk*a*  (58) 

Define  f^  as  follows. 

4  *(1-44)4  <69> 

[58,59]  4  =.  4+n  -  Vk  4  (60) 

2.6.2  Backward  Prediction 


Define  the  vector  composed  of  the  N+l  elements  leaving  the 
state  vector  at  the  kth  sample  time  (when  k/N  is  an  integer)  as  follows. 

T  A  r  .  1 

4  =  Luk-K  !-l  Uk-K  i-2  •*'  Uk-K  i-N  dk-D-NK2-N^  (61) 

The  linear  backward  prediction  of  pj<  based  upon  3-  is  Bk_N 
where  Bk  is  an  (N+l)  x  (K1+l+K2)  dimensional  matrix  of  backward  prediction 
coefficients.  The  backward  prediction  error  vector  is  defined  as 

h  =  -^k  "  Bk-N  ’  <62) 

and  the  backward  prediction  LS  cost  function  is 


f"  n-k/N  {  *  w  *  * 

backward  prediction  =  l  *  ^N^N^N*8^  * 

LS  cost  function  n=0 


Denote  the  optimal  value  of  this  cost  function  as  E°  and  define 
the  following  two  matrices. 


vb  t  k/jN  xn-k/N 


4iN  4iN  =  X  Vk-N  +  4  4 


,  ,  k/N 

Rub  =  I 
k  n=0 


.n-k/N 


4n4n 


=  X  R 


+  4  4 


The  derivation  of  the  backward  predictor  is  identical  to  that  of 
the  forward  predictor.  Because  the  definitions  given  by  Equations  (61) 
through  (65)  are  consistent  with  those  of  the  forward  predictor,  replace¬ 
ment  of  the  quantities  Fk,  f^,  v£,  e£,  r[,  j^+N  and  f^,  by  Bk,  b^,  Vk,  Ek, 
Rb,  ^  and  b^,  respectively,  into  Equations  (43)  through  (60)  produces  the 
mathematical  framework  of  the  backward  predictor.  The  expressions  needed 
for  subsequence  derivations  are  stated  below. 


\  B,.  =  \ 

k  k  k 


Bk  '  Bk-N  +  4  4 


2.6.3  Extended  State  DFE  Adaptation 


From  Equations  (39)  and  (61)  the  elements  of  and  pf<  are  never 
contained  within  the  same  state  vector.  To  establish  the  interrelation- 
ship  of  and  with  e£  and  vj ,  an  extended  DFE  is  considered.  The 
state  vector  of  this  extended  DFE  is  defined  below. 


-T 

x 

~k 


fu  u 
k  k- 


Uk-Ki-N  dk-D-N  dk-D-2N  ***  dk-D-NK2-N^ 


(70) 


Since  this  extended  DFE  is  just  the  original  DFE  with  the 
parameter  «i  increased  by  N  and  the  parameter  K2  increased  by  1,  Equations 
(8)  through  (38)  apply  here  directly.  Only  the  following  two  expressions 
related  to  the  extended  DFE  will  be  needed  however. 


[11] 


A  k/N 
=  1 
n=0 


n-k/N  -  -* 

-^nN  -SnN 


(71) 


\  3k  =  2k 


(72) 


where  gj<  is  the  extended  Kalman  gain  vector. 

2.6.4  Derivation  of  the  Algorithm 

To  relate  A|<  and  xj<  to  A|<  and  xj< ,  two  (K1+K2+N+2)  dimensional 
square  matrices  are  introduced. 

A 


It  is  easy  to  verify  that  Of  and  have  the  following  proper- 


(73) 


-  q;1 


0T  =  o'1 

%  % 


°f  4  = 


Qb  4  ■  - 


Next  the  extended  Kalman  gain  vector  is  expressed  in  terms  of 
the  forward  predictor  variables. 


°f  Sk  Qf  =  L  X"  k/N  Qf  ^iN  ^nN  Qf 


[75,77]  Qf  Ak  Q}  -  lQ  X 


n-k/N 


JjknN 


[11,44,  Q  A  qj  = 
48,78]  T  T 


Qf  \  i  =  H 


A  similar  derivation  is  now  performed  to  express  the  extended 
Kalman  gain  vector  in  terms  of  the  backward  predictor  variables. 

™  \  oj  ■  X  X"'MN  Qb  ^  Qb  (86) 

n=u 


T  n  t /N 

[86,76]  Q  S  QJ  ■  l  Xn'k/N  .  .  . 
b  k  b  n-0 

— nN 


[11,65,  Q  A  Qj  = 
66,87]  b  k  b 


*  *  | 
ViN  -nN  J 


Qb  \  4  =  Qb  4 


,,b* 

k  4 


[19,88]  QbAkQ  l 


[6^b8'  %\% 


k  k  k 


1  “■ 
■  v 


yy 


[97] 


(98) 


[62]  5,,  ■  £k  -  Bk_N  4 

[99]  Bk  =  (Bk.N  -  ^r1 
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[96] 

=  gk  +  Bk  uk 

[21] 

ek  =  dk-D  "  ^k-N  4 

[23] 

^k  =  ^k-N  +  -9k  Gk 

The  operation  of  a  T/N  fractionally-spaced 
Kalman  adaptation  algorithm  is  shown  in  Algorithm  2. 
that  the  direct  implementation  of  algorithms  such  as 

DFE  using  the  fast 
Experience  has  shown 

this  can  exhibit 

numerical  instability.  Often  adherence  to  the  IEEE  floating  point  stan¬ 
dard  results  in  a  stable  implementation.  Alternately,  the  algorithm  can 
be  restructured  in  the  form  of  the  so-called  square  root  filter  to  improve 
numerical  stability. 
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SECTION  3 

SIMULATION  DESCRIPTION 


This  section  describes  the  Monte  Carlo  computer  simulation  used 
to  evaluate  the  adaptive  equalizer  algorithms  discussed  in  thr  previous 
section.  The  simulation  is  partitioned  into  four  components  as  shown  in 
Figure  3.  This  simulation  structure  intentionally  resembles  the  signal 
flow  diagram  of  a  one-way  HF  link.  The  transmitted  waveform  is  first  syn¬ 
thesized  into  a  sampled-data  representation,  then  corrupted  by  the  propa¬ 
gation  medium,  and  finally  processed  within  the  receiver.  The  transmitted 
and  received  digital  data  streams  are  compared  to  accumulate  a  symbol 
error  rate  statistic.  Each  of  the  four  components  shown  in  Figure  3  are 
described  in  this  section. 

3.1  TRANSMITTED  SIGNAL  GENERATOR 

This  component  of  the  computer  simulation  generates  the  sample 
sequence  {d^  }  discussed  in  Section  2.1.  The  scope  of  this  report  is 
limited  to  binary  PSK  modulation,  so  all  channel  samples  are  binary. 
Further  the  information  bearing  symbols  (also  referred  to  as  bits)  are 
assumed  to  be  random,  independent  and  uncorrelated.  A  pseudo-random 
number  generator  implemented  in  software  generates  a  stream  of  random  ones 
and  zeroes  constituting  the  information  syr.iools. 

Recall  from  Section  ?  that  N  identical  and  consecutive  samples 
are  generated  for  pach  channel  symbol,  and  that  k  is  the  sample,  not 
symbol,  index.  Also  recall  that  for  BPSK  modulation  the  elements  of 


Figure  3.  Block  diagram  of  digital  PSK  HF  link  simulation. 


The  BPSK  receiver  that  does  employ  adaptive  equalization  is 
Receiver  Model  No.  2,  described  in  Section  3.3.2.  The  equalizer 
parameters  used  here  are  the  following: 


N  =  2 


j  1  x  io- 3  after  reception  of  an  information  symbol 
|  4  x  10-3  after  reception  of  a  training  symbol 

The  two  simulated  links  were  exercised  for  Es/N0  equal  to 
b  and  10  dB.  The  si gnal -to-noi se  ratio  for  the  DFE-assisted  modem  was 
actually  reduced  by  0.2948  dR  to  allow  for  the  seven  percent  training 
sequence  overhead.  Thus  the  two  modem  types  are  compared  on  the  basis  of 
equal  transmitted  energy  per  information  symbol. 


Where  E^  diminishes,  that  is  E^/E^  is  -°°  dB,  the  average 
bit  error  rate  for  Receiver  Model  No.  1  is  given  by 

P  =  1  Erfc  f/E  /N  ) 
e  2  ^  s  o; 

/here  Erfc(*)  is  the  familiar  complementary  error  function.  Evaluation 
of  this  expression  for  E^/N^  equaling  6  and  10  dB  yields  2.388  *  10-3 
and  3.8756  *  10-6,  respectively. 


The  average  bit  error  rate  for  Receiver  Model  No.  1  with  Es 
equaling  6  dB  is  plotted  versus  E^/E  in  Figure  9.  As  E^j  diminishes 
the  curve  appears  to  approach  the  predicted  value  of  2.388  x  10-  .  The 


RESULTS  AND  CONCLUSIONS 


This  section  presents  two  sets  of  results  that  were  generated 
during  this  investigation.  The  first  set  is  a  comparison  of  the  perform¬ 
ance  of  an  idealized  conventional  BPSK  receiver  to  that  of  a  DFE-assisted 
BPSK  receiver  in  the  context  of  the  time-invariant,  frequency-selective 
channel  model  described  in  Section  3.2.1.  These  results  are  presented  to 
illustrate  the  problem  of  intersymbol  interference  and  show  that  adaptive 
equalization  is  a  potential  mitigation  technique. 

The  second  set  of  results  show  the  DFE-assisted  modem  perform¬ 
ance  in  the  context  of  the  time-varying,  frequency-selective  channel  model 
discussed  in  Section  3.2.2.  These  results  give  an  indication  of  how  this 
modem  type  would  perform  over  a  range  of  doppler  and  delay  spreads 
exhibited  by  disturbed  HF  propagation  paths. 

4.1  PERFORMANCE  OF  RECEIVER  MODELS  NO.  1  AND  NO.  2 

IN  THE  CONTEXT  OF  CHANNEL  MODEL  NO.  1. 

The  Monte  Carlo  computer  simulation  was  exercised  to  investigate 
the  performance  of  the  two  BPSK  modems  described  in  Section  3.3  operating 
under  the  time-i nvari ant ,  frequency-selective  channel  conditions  discussed 
in  Section  3.2.1. 

The  BPSK  receiver  not  employing  adaptive  equalization  is 
Receiver  Model  No.  1,  described  in  Section  3.3.1.  It  does  not  use  a 
training  sequence  and  is  bit  synchronized  and  phase  locked  onto  the  signal 
energy  propagation  mode. 
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the  simulation.  Instead  the  decision  rule  offset  delay  parameter  D  is 
selected  so  that  all  signal  energy  associated  with  dk-D  is  in  the  feed 
forward  section  of  the  DFE  when  yk  is  generated. 

A  block  diagram  of  the  model  is  shown  in  Figure  8.  It  differs 
from  Receiver  Model  No.  1  in  three  ways:  1)  a  DFE  has  replaced  the 
N-sample  accumulation  function;  2)  the  parameter  D  is  larger  than  zero; 
and  3)  the  decision  symbol  dk_D  iS  fed  back  to  the  DFE  function. 

DFE  adaptation  can  be  implemented  using  either  Algorithm  1  of 
Section  2.3  or  Algorithm  2  of  Section  2.6,  with  one  modification.  The 
variable  d^_p  of  Step  4  of  Algorithm  1  and  Step  12  of  Algorithm  2  must 
be  replaced  by  d^.  The  decision  rule  for  dk  for  Receiver  Model 
No.  1  differs  from  Receiver  Model  No.  2  because  of  the  incorporation  of 
training  sequence  symbols  into  the  telemetry  format.  The  rule  is  given 
by: 

!+l  Real  {yk }  2.  0  and  ^k-D  was  generated  from  an 

information  symbol 

-1  Real  {y k }  <  0  and  dk-D  was  generated  from  an 

information  symbol 

dk-D  dk-D  was  generated  from  a  training  symbol 

3.4  ERROR  STATISTICS  ACCUMULATION 

This  part  of  the  simulation  simply  compares  the  sequence  {dk } 
with  {dk}  to  form  an  estimate  of  the  average  symbol  error  rate.'  The 
error  rate  estimate  is  calculated  by  dividing  the  number  of  incorrectly 
received  information  symbols  by  the  total  number  of  received  information 
symbols.  Training  sequence  symbols  are  only  used  to  adapt  the  DFE  and  are 
not  used  in  the  calculation  of  the  error  rate  estimate. 


3.3.1  Receiver  Model  No.  1:  An  Idealized  Conventional 

BPSK  Receiver  Model 

This  receiver  model  is  only  considered  in  the  context  of  the 
simply  delay  spread  channel  model  described  in  Section  3.2.1.  The 
receiver  is  assumed  to  be  synchronized  and  phase-locked  onto  the  mode 
propagating  signal  energy.  Since  the  signal  energy  mode  is  represented  by 
h°,  the  decision  variable  sequence  is  denoted  {d^}.  The  delay  parameter 
D  equals  zero  because  d^  can  be  generated  immediately  after  the 
kth  sampling  interval. 

This  receiver  model  is  represented  by  Figure  7.  A  value  for 
yi<  is  generated  (for  each  k  that  is  an  integer  multiple  of  N)  by 
simply  adding  u|<,  U|<_i,  ....  u^-n+i-  The  binary  decision  variable  d^  is 
then  formed  by  the  following  rule: 

1  Real  (yk}  >  0 

-1  Real  {y  }  <  0 

k 

where  again  d^  is  only  calculated  when  k  is  an  integer  multiple  of 
N.  By  definition,  the  channel  input  samples  d^ ,  dk_i,  ...,  d|<_N  +  1  are 
equal,  being  identically  generated  from  the  same  channel  input  symbol  as 
discussed  in  Section  3.1. 

3.3.2  Receiver  Model  No. 2:  A  DFE -Assisted  BPSK  Receiver 

This  receiver  model  is  used  for  both  channel  models  discussed  in  • 
Section  3.2.  The  DFE  explicitly  provides  the  phase  tracking  and  "fine" 
bit  synchronization  functions  that  are  assumed  to  be  perfect  in  Receiver 
Model  No.  1.  In  an  operational  receiver  a  "coarse"  bit  synchronization 
capability  would  be  required,  but  this  is  not  implemented  explicitly  in 
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Figure  6.  Delay  power  spectrum  of  Channel  Model  No.  2. 
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B  =  Doppler  bandwidth  (hertz) 

C  =  average  carrier  power  (watts) 

p  =  constant  selected  so  that  C  =  E  {h j^d ^ }  (note 
that  p  equals  1/2  for  N  less  than  four) 

For  the  purposes  of  this  report,  this  channel  model  will  be 
specified  by  the  four  parameters  C/Nq,  B,  N  and  T.  C  is  arbitrarily 
chosen  to  be  unity  without  loss  of  generality.  Since  the  parameter 
C/N0  has  units  of  watts/joule  or  hertz  rather  than  being  a  unitless 
ratio,  the  logarithmic  form  of  the  C/N0  specification  has  units  of  dB-Hz 
(decibels  relative  to  a  one  hertz  bandwidth).  The  delay  spectrum  for  the 
channel  model  is  illustrated  by  Figure  6. 

This  channel  model  shall  not  be  employed  to  investigate  conven¬ 
tional  BPSK  modems,  such  as  Receiver  Model  No.  1,  which  can  only  lock  onto 
one  of  the  four  modes.  Therefore  it  is  not  appropriate  to  designate  one 
mode  as  the  propagation  path  for  signal  energy  with  which  all  other  modes 
interfere.  Rather  the  received  energy  from  each  mode  of  Channel  Model 
No.  2  is  considered  signal  energy  even  though  the  four  modes  interfere 
destructively  with  one  another  on  the  average. 

3.3  RECEIVER  MODELS 

A  receiver  model  generates  the  decision  variable  sequence  |d|<_[)} 
(for  all  k  that  are  integer  multiples  of  N)  by  processing  the  channel 
model  outpu  .  symbol  sequence  {u|<}. 

Two  receiver  models  have  been  employed  for  this  investigation  - 
one  that  models  an  idealized  conventional  BPSK  receiver  and  another  that 
models  an  adaptive  OFF  assisted  receiver. 
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model  is  inadequate  if  the  modes  exhibit  different  doppler  frequency 
shifts.  Such  channels  are  said  to  have  a  doppler  spread  and  can  be  repre 
sented  by  a  channel  model  that  allows  a  distinct  phase  rotation  rate  for 
each  mode.  One  such  channel  model  is  described  below. 

3.3.2  Channel  Model  Ho.  2:  A  Simple  Delay  and  Doppler  Spread  Model 

This  model  can  also  be  represented  using  the  notation  of 

Section  2. 


L  =  13 

/Fp 

exp 

(0.375  • 

•  2ukBT/N } 

for 

m  =  0 

/c~P 

exp 

(0.125  ■ 

«  2-rrkBT/N } 

for 

m  =  3 

m  . 
hk  = 

/¥q 

exp 

(0.250  • 

•  2irkBT/N} 

for 

m  =  9 

/(Tp 

exp 

(0.625  < 

•  2iTkBT/N } 

for 

m  =  13 

^  0 

othe 

rwi  se 

E 


N  N/T 
o 


k  i 

n.  n .  = 

k  J  1 


k  «  j 

k  *  j 


where  k,  N,  n^,  N0  retain  the  same  definitions  as  in  Section  3.3.1 


N  =  ratio  of  sampling  rate  to  modulation  rate 

n^  =  complex  zero-mean,  additive  white  Gaussian  component 
of  Equation  2  of  Section  2.  Here  the  quadrature  real 
and  imaginary  components  of  n^  are  uncorrelated  and 
have  the  same  variance 

Nq  =  one-side  power  spectral  density  level  of  each  quadrature 
component  of  the  additive  white  noise  preceding  sampler. 

For  the  purposes  of  this  report,  this  channel  model  will  be 
specified  by  the  three  parameters  E  /N  ,  E  /E._r  and  N.  E  and  T  are 

S  0  S  1  o  1  S 

arbitrarily  chosen  to  be  unity  without  loss  of  generality. 

This  channel  model  can  be  visualized  using  Figure  5  which  shows 
a  two  mode  energy  delay  spectrum  resulting  from  the  transmission  of  an 
impulse.  The  earliest  (left  most)  spike  corresponds  to  the  mode  propa¬ 
gating  signal  energy  and  the  other  spike  corresponds  to  the  mode  propa¬ 
gating  intersymbol  interference  energy.  The  association  of  one  mode  with 
the  signal  and  the  other  with  the  intersymbol  interference  is  motivated  by 
how  this  channel  model  might  affect  a  conventional  PSK  receiver.  Assuming 
the  delay  between  the  two  modes  exceeds  a  channel  symbol  interval  (i.e.,  N 
is  less  than  20)  the  matched  filter  of  a  conventional  PSK  demodulor  cannot 
combine  energy  from  both  modes  that  arises  from  the  same  channel  symbol. 

This  channel  model  approximates  only  a  limited  range  of  HF  prop¬ 
agation  possibilities.  Three  examples  might  be  simultaneous  reception 
from:  1)  a  groundwave  and  a  skywave  propagation  path;  2)  a  single  hop  and 
a  double  hop  path;  and  3)  two  single  hop  paths  reflecting  off  distinct 
ionospheric  layers.  In  the  context  of  adaptive  equalization,  the  above 


effects.  Reference  16  reviews  the  current  state  of  the  art  in  HF  propaga¬ 
tion  analysis  in  the  context  of  nuclear  phenomenology. 


Since  the  current  understanding  of  HF  propagation  disturbances 
is  incomplete,  the  choice  of  a  channel  model  is  somewhat  arbitrary.  Here 
two  simple,  non-stati sti cal  channel  models  shall  be  employed,  which  have 
the  salient  characteristics  needed  for  this  investigation,  but  are  not 
general  representations  of  HF  channels.  The  first  is  time  invariant, 
exhibiting  delay  spread  but  not  doppler  spread.  The  second  exhibits  both 
delay  and  doppler  spreads. 

3.2.1  Channel  Model  No.  1:  A  Simple  Delay  Spread  Model 

Using  the  notation  of  Section  2,  this  channel  model  can  be 
described  by  the  following  set  of  parameter  values. 

m  =  0 

m  =  20  for  all  k 

otherwise 
k  =  J 
k  *  j 

E^  =  received  signal  energy  per  modulation  interval 

Ej^j  =  received  intersymbol  interference  energy  per 
modulation  interval 


L  =  20 


7T 


/T 


hk  ^ISl/T 


E  { n  n  * } 
k 


.1 


N  N/T 
0 


where 


k  =  sample  index 


{d|<}  can  be  taken  from  the  set  { 1 , - 1 } .  The  convention  used  in  the  simu¬ 
lation  is  as  follows:  the  binary  channel  symbol  0  is  mapped  into  N  con- 
sequetive  elements  of  { d^ }  having  the  value  +1;  and  channel  symbol  1  is 
similarly  mapped  into  N  consecutive  elements  having  the  value  -1. 

In  addition  to  the  correspondence  with  information  symbols,  some 
elements  of  {d^}  may  correspond  to  training  sequence  symbols.  A  train¬ 
ing  sequence  is  a  block  of  channel  symbols,  known  a  priori  to  the  re¬ 
ceiver,  which  are  periodically  injected  into  the  transmitted  data  stream. 
Training  sequence  symbols  are  used  within  the  receiver  to  initialize  and 
maintain  the  adaptation  algorithm  of  the  equalizer,  as  will  be  discussed 
in  Section  3.3. 

The  telemetry  format  of  the  simulated  HF  link  is  shown  in 
Figure  4a.  The  symbol  stream  is  partitioned  into  telemetry  blocks  of 
l_TB  symbols  in  duration.  Each  telemetry  block  starts  with  the  same 
sequence  of  Lj$  training  symbols  (represented  by  cross-hatching  in  the 
figure),  followed  by  a  Ljb-Lts  length  segment  from  the  information 
symbol  stream.  The  multiplexing  of  the  training  sequences  into  the  data 
stream  is  represented  by  Figure  4b.  The  generated  sequence  {d|<}  is 
taken  as  input  by  the  channel  models. 

3.2  CHANNEL  MODELS 

Much  experimental  and  analytic  work  has  been  done  to  charac¬ 
terize  HF  channels  (e.g..  References  7-15),  but  universally  accepted 
models  of  HF  propagation  disturbances  have  yet  to  be  developed.  In 
natural  environments,  a  wide  range  of  channel  cha racteri sti cs ,  such  as 
absorption,  delay  spread  and  doppler  spread,  have  been  observed  and  the 
most  severe  effects  tend  to  occur  on  auroral  propagation  paths.  An  impor¬ 
tant  class  of  artificially  induced  effects  are  those  resulting  from  high 
altitude  nuclear  bursts.  Such  effects  are  not  presently  well  understood, 
but  are  estimated  to  be  at  least  as  severe  as  naturally  occurring  auroral 
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bit  error  rate  increases  with  the  level  of  intersymbol  interfere  as  would 
be  expected.  The  severity  with  which  this  error  rate  increases  accounts 
for  the  rarity  with  which  this  modem  type  is  employed  in  HF  communication 
links. 


Two  approaches  for  improvi;  g  the  average  bit  error  rate  were 
explored:  increasing  the  signal -to-noi se  ratio;  and  substituting  Receiver 
Model  No.  2  for  Receiver  Model  No.  1.  Figure  10  presents  a  comparison  of 
the  results  of  these  approaches  with  the  results  of  Figure  9. 

Consider  first  the  difference  in  error  rate  characteristics  when 
Es/N0  is  increased  from  6  dB  to  10  dB.  As  previously  mentioned,  this 
would  result  in  a  decrease  in  average  bit  error  rate  from  2.388  *  10-3  to 
3.8756  x  10-6  for  E j 5 £  equal  to  zero.  Figure  10  shows,  though,  that  as 
E .  , /E  is  increased,  the  10  dB  error  rate  characteristic  rapidly 
approaches  the  same  degraded  performance  level  experienced  for  6  dB. 
Apparently  as  E j 5 1 /E s  is  increased,  the  significance  of  the  additive  white 
Gaussian  noise  on  the  error  mechanism  diminishes  in  comparison  to  the 
significance  of  intersymbol  interference.  Therefore  increasing  the 
transmitted  carrier  power,  which  increases  both.E^j  and  F^  but  leaves 
E 1  s I /E s  una^terecl»  does  not  effectively  combat  the  effects  of  high  levels 
of  intersyrnbol  interference. 

Now  consider  the  difference  in  error  rates  with  Es/N0  held 
at  6  dB  when  Receiver  Model  No.  1  is  replaced  with  Receiver  Model  No.  2 
and  the  attendant  change  in  telemetry  format  to  accommodate  a  training 
sequence.  When  Ej^j/E  is  very  small.  Receiver  No.  2  suffers  an  effective 
0.4  dB  loss  relative  to  the  asymptotic  error  rate  exhibited  by  Receiver 
No.  1,  but  as  E j ^ j /E ^  increases  the  error  rate  for  Receiver  No.  2 
decreases,  rather  than  increases.  The  0.4  dB  loss  in  the  absence  of 
intersyrnbol  interference  is  due  in  part  to  the  0.2948  dB  training  sequence 
overhead  mentioned  earlier,  and  in  part  to  the  random  misadjustment  of  the 


Bit  Error  Rate 


Ratio  of  E,s,  to  Es  (dB) 


DFE  coefficients  caused  by  the  influence  of  additive  white  Gaussian  noise 
on  the  adaptation  algorithm.  The  decrease  in  average  bit  error  rate  with 
increasing  ^ j s I s  resu^s  ^rom  the  DFE  using  the  intersymbol  interference 
energy  as  additional  signal  energy. 


4.2  PERFORMANCE  OF  RECEIVER  MODEL  NO.  2  IN  THE  CONTEXT 

OF  CHANNEL  MODEL  NO.  2. 

The  Monte  Carlo  simulation  was  exercised  to  investigate  the 
performance  of  the  DFE-assisted  BPSK  modem  described  in  Section  3.3.2 
under  the  time  and  frequency  selective  channel  conditions  discussed  in 
Section  3.2.1.  Channel  Model  No.  1  is  time-invariant,  allowing  the 
presentation  of  the  results  in  Section  4.1  to  be  in  a  form  that  is 
independent  of  the  data  rate. 

Here  Channel  Model  No.  2  is  time-varying,  and  it  is  convenient 
to  consider  one  specific  information  symbol  rate,  which  was  arbitrarily 
chosen  to  be  2400  baud.  The  DFE  parameters  used  to  generate  all  the 
results  presented  in  the  section  are  as  follows: 


LTB  "  96 
N  =  2 

«!  =  8 
K2  =  16 


|2  x  10- 3  after  reception  of  an  information  symbol 
| 8  x  10- 3  after  reception  of  a  training  symbol 


T  =  1/4800  secopds 


Receiver  Model  No.  1  is  not  considered  here  because  it  does  not 
have  a  phase  tracking  capability.  The  investigation  of  Receiver  Model 
No.  1  in  the  context  of  any  time-varying  channel  model  would  be  pointless. 

The  stochastic  gradient  algorithm  discussed  in  Section  2.3  was 
used  within  Receiver  Model  No.  2  to  adapt  the  DFE  coefficients.  This 
algorithm  has  an  inherent  phase  tracking  capability  which  is  similar  to 
that  of  a  first  order  phase  lock  loop.  Since  each  mode  has  a  fixed 
doppler  offset,  the  adaptation  algorithm  must  track  a  phase  ramp.  This 
results  in  a  lag  error  for  each  DFE  coefficient,  and  the  adaptation 
algorithm  could  have  been  redesigned  as  a  second  order  phase  tracking  loop 
to  reduce  these  lag  errors.  However,  by  not  specifically  optimizing  the 
receiver  for  Channel  Model  No.  2,  the  results  obtained  here  are  indicative 
of  the  performance  of  this  modem  type  for  other  channels  of  comparable 
doppler  spread  but  exhibiting  more  random  phase  fluctuations. 

The  results  take  the  form  of  average  bit  error  rates  plotted 
versus  either  the  average  carrier  power-to-noi se  spectral  density  level 
C/N0  of  the  doppler  bandwidth  8. 


Figure  11  shows  the  bit  error  rate  versus  C/N0  for  B  fixed  at 
four  hertz.  As  expected,  the  error  rate  decreases  monotoni cal ly  with 
increasing  C/N0.  However,  the  error  rate  characteristic  flattens  for 
values  of  C/N0  above  46  dB-Hz.  When  C/NQ  is  below  this  value  most 
detected  bit  errors  are  caused  by  the  AWGN  component  of  the  channel  model 
output.  As  C/N0  is  increased  the  AWGN  component  diminishes  in  compari¬ 
son  to  the  residual  intersymbol  interference  which  the  adaptive-DFE  does 
not  correct.  Further  increasing  C/N0  does  not  appreciably  improve  the 
error  rate  because  the  effects  of  intersymbol  interference  are  not 
reduced. 


error  rate  versus  C/N0  for  Recel 
ler  bandwidth  of  4  hertz. 


Figure  12  shows  the  error  rate  versus  C/N0  characteristic  for 
several  coppler  bandwidth  values.  Figure  13  shows  the  same  data  plotted 
in  the  form  of  error  rate  versus  B  characteristic  curves  for  several 
values  of  C/N0.  The  DFE  coefficients  misadjustment  due  to  adaptation 
lag  increases  with  B,  and  the  resulting  effects  upon  the  error  rates  are 
shown  in  both  of  these  figures. 

4.3  CONCLUSIONS 

The  adaptive  DFE  is  an  effective  solution  to  intersymbol  inter¬ 
ference  if  the  channel  impulse  response  is  constant  or  very  slowly 
changing.  For  severe  levels  of  intersymbol  interference,  the  inclusion  of 
adaptive  equalization  into  the  link  design  was  shown  in  Section  4  to  be 
much  more  beneficial  than  increasing  the  received  signal -to-noise  ratio. 

Th<  effectiveness  of  the  adaptive  DFE  can  be  limited  by  the 
rapidity  of  channel  fluctuations,  parameterized  in  Section  4.2  by  the 
doppler  bandwidth  B.  The  specific  2400  baud  modem  considered  in  that 
section  exhibited  usable  channel  capacity  for  values  of  B  up  to  4  Hz.  For 
C/N0  equal  to  44  dB-Hz  and  B  equal  to  4  Hz,  the  measured  channel  symbol 
error  rate  was  about  0.04.  To  make  such  a  modem  useful  for  practical 
applications  one  would  employ  some  form  of  channel  symbol  error  control. 
Augmenting  this  modem  with  rate  1/2,  constraint  length  7  convolutional 
coding,  interleaving  and  Viterbi  decoding  would  produce  a  1200  band  data 
link  with  an  expected  error  rate  of  less  than  10"3.  Employing  rate  1/4  or 
rate  1/8  codes  instead  of  the  rate  1/2  code  would  yield  600  and  300  baud 
links  with  much  lower  error  rates. 
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Figure  12.  Bit  error  rate  versus  C/N0  for  Receiver  Model  No.  2  for 
various  doppler  bandwldths. 


Figure  13.  Bit  error  rate  versus  doppler  bandwidth  for  various 
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